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Abstract 

The interlayer magnetoresistance of a quasi-two-dimensional layered metal with a J-wave pseudogap is 
calculated semiclassically. An expression for the interlayer resistivity as a function of the strength and direc- 
tion of the magnetic field, the magnitude of the pseudogap, temperature, and scattering rate is obtained. We 
find that the pseudogap, by introducing low-energy nodal quasiparticle contours, smooths the dependence 
on field direction in a manner characteristic of its anisotropy. We thus propose that interlayer resistance 
measurements under a strong field of variable orientation can be used to fully characterize an anisotropic 
pseudogap. The general result is applied to the case of a magnetic field parallel to the conducting layers 
using a model band structure appropriate for overdoped T(22Q\. 

PACS numbers: 71.27.-Ha,72.10.Di,71.10.Ay 
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INTRODUCTION 



High temperature superconducting cuprates, organic charge transfer salts, some heavy fermion 
materials and a host of other intriguing electronic systems, are layered metals in which electrons 
are approximately confined to a given atomic layer. Much of the interesting behavior of these 
materials arise because of strong electronic correlations within a single layer. Surprisingly, it 
turns out that one of the most effective means of accessing in-layer properties, particularly those 
properties that are highly anisotropic within a layer, is to measure interlayer electronic transport 
coefficients in a strong magnetic field[[l[y,[3[|4, 5, 6, 7, 8]. 

The interlayer electrical resistivity p,, depends on the direction of the magnetic field in a man- 
ner that is highly sensitive to the anisotropy of the quasi-2D band structure. High-resolution 
maps of the Fermi surface, and other band structure properties, have already been obtained by 
fitting data to calculations based on semiclassical magnetotransport theory. This technique 



has been applied to a wide variety of layered materials including overdoped cuprates I 
ruthenatesDSl l6|], and organic charge transfer salts itzl Isl]. The p^^ data also contains information 
about in-plane scattering and can be used to study the directional dependence of elastic and in- 

11511 . Notably, it has been used to reveal a T-linear, 



elastic scattering rates |19L [10, 1 111 . 



anisotropic scattering contribution in overdoped cuprate superconductors that appears to be tied to 



superconductivity itself ISllailTD. It is important to press further, to ask what other anisotropic 



properties of the metallic layers can be detected and characterized via interlayer transport in high 
magnetic fields. 

In this article we ask what interplane transport data can tell us about an anisotropic pseudo- 
gap Ajj in quasi-2D metals. Since an anisotropic gap in the density of states will aff"ect the field- 
direction dependence of p,^, we expect that interlayer magnetoresistance can be used to map out Ak 
as well. A natural application of this technique would be to slightly overdoped cuprates. For these 
materials, a model of a 2D metal with a small J-wave pseudogap (that is starting to emerge with 
reduced doping) is a plausible description of the metallic state at fields above Hc2 and semiclassi- 
cal calculations of p,^ may adequately capture transport properties. To extract from p,^ information 



about the doping, temperature and field dependence of A^ would be of great va 



standing the relationship between the pseudogap and superconductivity [| 1 8 



ue towards under- 



19h . The effects of a 



non-zero At may already be present in existing interlayer resistance data on slightly overdoped 
cuprates, convoluted with the effects of anisotropic scattering [20]. If so, a reinterpretation of these 
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data using models that incorporate a pseudogap could be fruitful. 

We study a model with well-defined electronic quasiparticles existing in the presence of a d- 
wave pseudogap in the density of states. The manner in which the opening of the pseudogap will 
change the interlayer resistivity is predicted and the following main results obtained: 

i) An expression for the interlayer resistance p„ in the semiclassical limit in a strong magnetic 
field of arbitrary strength and direction. 

ii) For the simple case of a field parallel to the layer, with arbitrary intralayer orientation (pB, the 
quantitative effect of a pseudogap on p^zispB) is calculated using a realistic model band structure. 
The average magnitude of Pzz{(Pb) varies non-monotonically with the size of the pseudogap while 
its (pB dependence is modified in a manner distinctive of the pseudogap symmetry. A strongly 
anisotropic normal-state p^zifp) is smoothed by the pseudogap through the introduction of new 
low-energy current contributions associated with J-wave nodes. 

Considering our results in light of the success of the AMRO technique in extracting band struc- 
ture and scattering parameters of cuprates, we propose that this technique should also prove to be 
a viable means of obtaining a T-, B- and doping-dependent parametrization of the t^-wave pseu- 
dogap. 



SEMICLASSICAL PICTURE OF PSEUDOGAP STATE 



As a simple model of the J-wave pseudogap state one can use the normal (diagonal) part of the 
BCS Green's function, taking the anomalous part equal to zero. The Green's function is 

2 2 

Go(a>, k, X) = + -\-] (1) 

\a> - Ek (x> + E\^l 

with band energy pseudogap Ak and relative spectral weights for the electron and hole terms: 

"k = ^(l+^k/£k) , v2 = i(l-^k/£k) (2) 
and a quasiparticle energy given by 

£k = + (3) 



Recently Yang, Rice and Zhang (YRZ^ HllN proposed an ansatz for the coherent part of the 
Green's function in the pseudogap state of high-temperature superconducting cuprates based on 
renormalized mean field theory calculations of the resonating valence bond state in the ? - 7 model. 



It has a similar form to Eq. [H differing only by the appearance of small additional terms in the band 
energy that break the particle-hole symmetry (also, YRZ proposed specific doping dependencies 
of the overall magnitude of the spectral weight of the coherent part, as well as of the band hopping 
parameters and Aq). YRZ found good agreement between their model and several characteristic 
features of ARPES data[2^- The interlayer resistance of the YRZ model could be studied using 
the same approach followed in this article. We consider here the more familiar BCS expression 
(Eq. [T]) in order to illustrate the qualitative changes to the normal-metal p,^ that are induced by 
turning on A^. 

Eq. [T]can be viewed as a description of a two-band metal with band energies of i^k, measured 
from the Fermi level, and k-dependent spectral weights. At T = the lower band is filled, the 
upper band empty and their nodal crossing point lies exactly at the chemical potential. If the 
imaginary part of the self-energy correction to Eq. [T] is small (compared to relevant co) then the 
quasiparticles in each band are well-defined and transport properties can be calculated using a 
semiclassical Boltzmann approach. 

For the semiclassical picture to be applicable, the quasiparticles in each band must remain well- 
defined, i.e. the imaginary parts of the self energy correction to Eq. [Umust be small compared to 
relevant frequencies oj. At low temperature and frequency, impurity scattering will dominate. The 



associated scattering rate can be obtained following the procedure for d-wave superconductors [1231] 
and it is known that, at sufficiently low frequency the impurity scattering rate becomes larger than 
the frequency so the semiclassical picture of transport is not applicable. At high temperature and 
frequency, strong inelastic scattering will also render the semiclassical approach invalid. However 
there may exist an intermediate frequency range for which both the impurity and inelastic scatter- 
ing rates are relatively small. In this range, quasiparticles are sharply defined and the scattering 
rate r'^co, k) can be evaluated at the quasiparticle pole co = Ek- We assume that such a frequency 
range exists and calculate the interlayer resistivity in a magnetic field using Boltzmann theory. 

INTERLAYER RESISTANCE IN THE PSEUDOGAP STATE IN THE PRESENCE OF AN ARBI- 
TRARY MAGNETIC FIELD 

To have interlayer current there must be a finite amplitude t± for hopping between adjacent 



layers. However, according to Kennett and McKenzie [|lO|] . the form of the interlayer conduc 



tivity does not depend on whether or not interlayer transport is coherent (i.e. it does not de- 
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pend on the relative magnitude of and fiT'^) as long as in-plane momentum is conserved 
during interlayer hopping. We may thus carry out the calculation of the interlayer conductiv- 
ity by supposing that a 3D Fermi (quasi-cylindrical) Fermi surface exists even when we are in 
the regime in which the Bloch vector in the interlayer direction is not well-defined. Taking 
advantage of this, we simply add to a term -2tj^{k^,ky)co^{k^c) where tj_{k^,ky) is the inter- 
layer hopping coefficient and c the distance between layers. The associated interlayer velocity is 
v^(kx,ky,k^) = 2cfi~^tj^(kx,ky) sin(fc^c). The calculation of the interlayer current is done to lowest 
order in v^. 

The Boltzmann equation in a weak electric field 3 along the z axis and a magnetic field B of 
arbitrary strength and direction is: 



where the total distribution is / = /o -I- /i with fi = -{dfQldE^)g, fo{x) is the Fermi function 
and I[g]js the collision functional. The auxiliary time variable t is defined by the equation of 



motion|124|] 



— = -e\g X B (5) 



dk 

dt 

where = dE^/dk. Eqs. [5] and H] are solved to obtain the distribution function, which is inserted 
into the expression for the interlayer current: 



Jz(t) = ^ J dkvMt])MtMt]). (6) 

The current is found by taking a ^-Fourier transform of j^{t) and evaluating in the zero-frequency 
limit. The spectral weights from the two bands combine simply to give + = 1 for the particle- 
hole symmetric case. 

For a field B = 5(sin 6^ cos (ps, sin 6^ sin 0, cos Ob), Eq. [5]gives dcpjdt = codE, (p, 9b) where the 
cyclotron frequency is 

HcociE, (f>, Ob) = eB cos dB^^^—^. (7) 

The cylindrical (p variable parameterizes the cyclotron orbit around a closed energy contour E^^ = 
E. Any point on the projection of this contour onto the k^ - ky plane is written as = 
ksicpXcos (f>, sin (p) where kE{(p) is the radial cylindrical distance measured from some arbitrary 
point in the region enclosed by the contour. In the normal state we can use a single energy contour 
(the Fermi surface) k^ = ky = kf{(p){cos (p, ^m(p). 
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The momentum varies according to 

= - tan OBj^\kE{(p) cos(0 - 0^) j, (8) 

which results in a periodic oscillation of the interlayer velocity v^{k^\t\) that is determined by the 
direction of the field angle Ob- 

Finally, since is independent of k^ in the collision functional, the integral over h of the 
'scattering in' term vanishes by symmetry (this is true not only for scattering from point de- 



fects. 



DUt 



'or any other scattering mechanism that can be regarded as spatially confined to a single 



plane [|l2ill5|]). We are left with a relaxation time-description: I\_g^ = -gk/T(Ek) with the current 
relaxation rate equal to the total quasiparticle scattering rate T~\a) = Ek). The fact that vertex cor- 
rections vanish to lowest order in in the calculation of the interlayer resistivity is a considerable 
simplification. It means that we can use any appropriate model for the scattering rate, including 
elastic or inelastic scattering or even a sum over several different mechanisms. 

We insert these expressions into Eq. IH formally solve for gk(t) and use this in Eq. [6] to obtain 



the interlayer conductivity cr^^ = l/p^ 



f 

%J — o 



^^^(B)=A| '^Jt^I^ (9) 



^ r^Ml r il^Mi,M (10) 

Jo ^OJc{4>2) fiojcicpi) 



where 



A = —rceB cos 9b, 



M(0i,(^2) = G((^i,(^2)COSO[0i,02] 

d(f>' 
coc{(p')r{(f>) 



G(0i,<^2) = exp(- —X—-) 

J<A, 



O(0i,^2) 

— — = fe(0l) COS(0i - 0b) - kE(cf>2) COS(02 - ^b)], 

c tan 6b 

and P = G(0, 2n). The function G(0i, ^2) is the probability that a quasiparticle can proceed from 
01 to 02 along its cyclotron orbit without being scattered so P is the probability that a quasiparticle 
completes an orbit. 

Eq. |9]has been written in the same form as the corresponding expression for a normal metal |lo| . 
However a»c(0), r, ^^(0) and P all depend on energy in the pseudogap state (though we have not 
always written this explicitly). Moreover, the and ^^(0) variables must be interpreted differently 
in this expression depending on whether the energy is greater or less than max A^. This is because 



these variables are defined with reference to a closed 2D cyclotron orbit but the orbits (i.e. the 
energy contours) have different topologies depending on the relative size of E and max Ak as shown 
in Fig. [U For E < max Ak there are four equivalent (banana- shaped) contours closed around nodal 
points so a node can be taken as an orbit center with the polar angle (p parameterizing position 
along the contour. Thus kE{(p), which is measured from the node to the contour, depends strongly 
on both E and (p. (We should include an overall sum over the four nodes in Eq. |9l though this has 
not been written explicitly. There is no mixing of different nodes since an electron remains on a 
single nodal contour during cyclotron motion and the current contribution from each nodal region 
can be obtained separately.) For E > max Ak a single contour encircles the entire normal state 
Fermi surface and kE{(p), measured from a central point, is weakly anisotropic-i.e. its anisotropy 
is that of the normal state Fermi surface. The energy integral in Eq. |9]must be broken up into low 
and high energy regions with the ^^(0) variable defined accordingly. 

Eq. |9] is the main result of this article. This expression could be used in fitting procedures 
similar to those applied in the normal state of overdoped cuprates. The magnitude of the pseu- 
dogap as, say, a function of doping, temperature and field strength in overdoped systems could 
then be extracted. A typical set of fitting parameters might include hopping amplitudes describing 
the normal state band structure and interplane hopping magnitude ?x (the values of which would 
be constrained by independent measurements and would be expected to be independent of tem- 
perature and weakly dependent on doping), the normal-state scattering rate r"' (which can also 
be independently estimated) and the gap magnitude. Additional parameters could be incorporated 
if one were to go beyond the nearest neighbor expression for Ak, or to include anisotropy in the 
scattering rate. Overall, the number of parameters would not have to exceed that used in previous 
normal state analysis. 

We will not undertake detailed numerical evaluations of Eq. |9]in this article but will discuss, 
in the remainder of this section, some of the general features of this expression that distinguish it 
from the familiar normal state result. The contribution to the conductivity, Eq. |9l that comes from 
energies E » max Ak will be identical to the normal state expression. So, the total conductivity 
is a weighted sum of the normal state value and the low-energy (i.e. ^k < max Ak) contribution 
associated with the pseudogap. The relative weighting is controlled by the value of Ao/^bT. The 
properties of the low-energy (pseudogap) contribution to the conductivity are qualitatively dif- 
ferent than those of the high-energy (normal state) contribution. It is more strongly temperature 
dependent and is less sensitive to the direction of the magnetic field. 
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To get some feel for the low-energy contribution to the magnetoconductivity, associated with 
the pseudogap, we consider the nodal limit E « Aq for which Ar and can be linearly expanded 
about nodal points. In the nodal limit we write A^ = V2k2 and = Vfki where ^2 and ki are 
momenta parallel and perpendicular to the Fermi surface, respectively. The radius of the energy 
contour with energy E is given by 

hkEiE,(f>) = — ^ (11) 

yjvj cos2(0 - 0„) + vl sin^(0 - (pn) 



and the cyclotron frequency is 



a)c(4>,E) = eB cos 6bE ^ 



vl cos^(0 - 0„) + vj sin^(0 - 0„) 



(12) 



where 0„ is the direction of the node (0„ = +n/4, ±3n/A). Since Vj >> the energy contour is a 
narrow ellipse and the cyclotron motion of the quasiparticle slows down dramatically as it crosses 
the Fermi surface ki = 0. 

If the probability P is small, then quasiparticles are unlikely to complete cyclotron orbits with- 
out being scattered and the field-dependence of the conductivity is weak. The field-dependent 
effects of interest (i.e. the sensitivity of cr^, to in-layer anisotropy and the AMRO) occur when P is 
of order 1 . The quantity P depends on the scattering mechanism and, generally in the pseudogap 



state, on energy E. In the simple case of point defects, the scattering rate ll23ll is approximately 
given by: 

T-\E) = T-o'[v(E)/voV (13) 

where ' is the normal state scattering rate, v{E) and vo are the densities of states in the pseudo- 
gap and normal states, respectively, and rj = +1 (or -1) in the Bom (or unitary) limit. For unitary 
scattering (to which we henceforth restrict ourselves) there is a cancelation in factors of the quasi- 
particle density of states so that P becomes energy-independent. In this case, P has roughly the 
same value as in the normal state. So, in strong fields, we can ignore the effects of scattering (i.e. 
set G{(f)\,(p2) = 1) in both the high-energy (normal state) contribution and the low-energy (nodal 
limit) contribution to a,^. This simplifies the following discussion. 

Both the sensitivity of cr^^ to the anisotropy of the 2D band structure and the AMRO effect 
originate from the argument O(0i, ^2) of the cosine in Eq. [9l The cosine oscillates rapidly when 
kfC tan 6b is large, and kills the integral everywhere except at special momentum directions, which 
depend on field orientation 0^. As discussed in Ref. [l^ the conductivity is thus dominated by the 
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small region where both 0i and (p2 are close to a special direction defined by the solution of 



d(p _ 



kf{4)) cos(0 - ^b) 



0. (14) 



Since the field direction determines the value of (pi and 02 that dominate the integrals, the band 
structure parameters are evaluated at a symmetry-unique point on the Fermi surface that can be 
tuned by field direction, allowing the Fermi surface to be mapped out. Also, since the scale of the 
rapid oscillation is set by kfC tan Og, the overall magnitude of the conductivity oscillates in 6b with 
a period determined by this quantity (this is AMRO). 

However, when we apply this reasoning to the low-energy pseudogap contribution, we find that 
such strong dependence on field direction angle is not expected. The solution to Eq. [14] in the 
nodal limit is - 0„ = arctan[Q'^(0B - 0„)] where a = v//v2 >> 1. The large factor means that 
the dominant value of - 0„ will almost always be close to n/l, i.e. close to the point at which 
the nodal energy contour crosses the Fermi surface, independent of the direction of field. (The 
only exception would be if the magnetic field were pointed precisely in a nodal direction.) So, the 
dependence on the field direction (ps is far weaker in the low-energy pseudogap contribution than 
it is in the normal state. Moreover, the scale for the oscillatory dependence in the nodal limit is 
fc^ctan^B « (£'/v2)ctan0. For energies E « Aq this quantity will be much smaller than one for 
any 6b + n/l. The argument of the cosine in Eq. |9]will be small and no oscillatory dependence on 
field angle 6b will be seen. Even at temperatures as high as kBT /Aq « 1 we do not expect to see 
prominent AMRO coming from the pseudogap contribution to the conductivity. This is because 
the integral over energy will average fc^ over all values from to nearly kf, giving no sharp period 
for oscillatory behavior. 

These qualitative arguments suggest that the low-energy (pseudogap) contribution to the con- 
ductivity will not show the strong field-direction dependence characteristic of the high-energy 
(normal state) contribution. (A detailed analysis is needed, however, to account for the strong 
energy-dependence of the scattering rate that could change this picture by giving dominant weight 
in the integral to a particular energy range.) So, the onset of the pseudogap should have the generic 
efi'ect of smoothing the dependence on field angle. Nevertheless, this smoothing will proceed in 
a particular manner that is characteristic of the anisotropy of the pseudogap. None of the above 
effects would occur for an isotropic pseudogap, and the J-wave case discussed here could be dis- 
tinguished from alternative forms since the arrangement of nodal points would have a different 
relationship with the normal state band anisotropy. 



In the next section we consider the simple limit of a field in the layers, i.e. 6b = n/2. This 
is done to provide a more quantitative description of the effect that a pseudogap has on the field- 
direction anisotropy of p^{^b)- Also, theoretical expressions for p^, with which we can compare 
our results, have been obtained previously using a different formalism. 



CASE OF A FIELD PARALLEL TO THE LAYERS 

The general result Eq. |9]can be evaluated in the limit 63 — > n/2 (i.e., for the case of a field 
in the layers) by employing a stationary phase approximation but it is simpler to go back to the 
beginning of the derivation and make this assumption. When B is in the layers: 

where k is the momentum in the plane and 

Qc(k) = ec\Vg X B|. (16) 

In the normal state, t'\E) = t'^ and Qc(k) = are both independent of energy so the integral 

over Ek gives unity and Eq. [15] reduces to a Fermi surface average. The magnetic field becomes 
important when the 0-averaged quantity Q.c ~ ecv/B becomes comparable to the scattering rate 
1 /r. Note that the criterion for field effects Q.ct ^ 1 is more favorable by a factor of kfC than the 
corresponding criterion for in-layer transport, where kfC « 10 is typical in cuprates[4|]. 

Eq. \T5\ can also be obtained using a tunneling Hamiltonian approach, and a similar result 



was thus obtained in Ref. |29|. The tunneling current is expressed as a convolution of spectral 
functions on adjacent layers. The gauge can be chosen such that the difference in the vector 
potential between adjacent layers is A = c{By, -B^, 0) and the corresponding spectral functions 
differ only by a momentum shift equal to eA. Evaluating the spectral functions in the quasiparticle 
approximation and using Clc * (£^k-eA-£^k) one obtains Eq. \T5\ The advantage of the semiclassical 
approach followed here is that it can be generalized to describe fields out of the layers (Eq. |9l). In 
the remainder of this article we will, however, focus on the simple case of Eq. \T5\ We go beyond 
the ksT « Aq nodal limit considered in Ref. 29 to consider arbitrary UbT I Aq and a realistic 
normal-state band structure for cuprates in order to study the effect of a small pseudogap on the 
(pB dependence of p,^. 

In a strong magnetic field, Qct >> 1 so Q.c{4>b)t » 1 at typical (pB, the Fermi surface average 
in Eq. [T5]is dominated by the k values for which Q.ci^) = 0> i-e. by k for which the quasiparticle 
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FIG. 1: A small pseudogap reduces the dependence of p^^ on the direction of a magnetic field B parallel 
to the layers. Upper inset: The dashed (green) curve is a Fermi surface, closed around the comer M point 
of the square Brillouin zone, and the hatched curve indicates the magnitude of the <i-wave pseudogap. Main 
panel: When the field is large, the interlayer current is dominated by k-points on low-lying energy 
contours at which the electron velocity = dEk/dk is parallel to B. The solid curves show low-lying 
^k-contours (moving outward from the node, the contours are for E/Aq = 0.05,0.5, 1, 1.5,2,2.5,3) in the 
upper right quadrant of the M-centred Brillouin zone. The arrows (each parallel to B) are located at the 
dominant k-point for each contour. In the normal state, the dominant k is the point where the largest (green) 
arrow intersects the dashed (Green) Fermi surface. In the pseudogap state, the dominant k are spread over 
a large range that extends from the normal state point (for E^ > max Ak) to the node (for E\i « max Ak). 
The opening of a pseudogap effectively spreads the current contribution over the Fermi surface, thereby 
smearing the 0b dependence of p^^. 

velocity is parallel to B. This means that the normal-state interlayer resistivity is determined by the 
values of band parameters at a particular point on the Fermi surface k = k* = kf {(()*)(cos (p* , sin 0*) 
where the value of 0* is controlled by (ps (in an isotropic system 0^ = 0*). Moreover, the resistivity 
is independent of t ' in strong fields since the current is limited by classical magnetoresistance 
rather than scattering. Upon varying (ps, one can use p„ to effectively map out the ^-dependence 
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of the in-plane band parameters. 

In the pseudogap state, the energy dependence of flc(k) changes this simple picture, as illus- 
trated in Fig. [TJ For energies » maxA^, the energy contours of the pseudogap state are 
almost identical to the Fermi surface itself. So the contribution to p^^ that comes from energies 
much larger than max are the same as in the normal state. However, when a pseudogap opens 
up (i.e. once Aq becomes comparable to ksT) the conductivity begins to receive significant con- 
tributions from energies < max Ar. The associated low-energy energy contours are centered 
on nodes and the k point on such an energy contour where Cld^) vanishes is far removed from 
the corresponding normal state point k*. This means that a small pseudogap results in contribu- 
tions to p,, coming from a much broader range on the Fermi surface, thereby weakening the 0g 
dependence. 

This loss of 0B-dependence occurs initially without a corresponding increase in the magnitude 
of cTj.^. In fact, the effect of turning on a small pseudogap in the presence of a strong in-layer 
magnetic field is to decrease the (pB averaged interlayer resistivity, as shown in Fig. [21 This comes 
about because classical magnetoresistance is relieved by the pseudogap, both through a reduction 
of the average quasiparticle velocity and through the increased range of k that contribute to the 
current. The effect is independent of r ' at large fields, as noted above, so the energy dependence 
of the scattering rate (i.e., whether we are in the Born or unitary limit) does not matter. For 
a sufficiently large pseudogap, the reduction of the carrier density overcomes this effect, so 
reaches a minimum at Ao/kBT x 1 and thereafter increases, eventually becoming very large for 
Ao/ksT » 1 when the current comes only from the nodal regions. 

CALCULATION OF INTERLAYER RESISTIVITY USING MODEL BAND STRUCTURE OF 
T^2201 

To obtain a more quantitative picture of the (pB dependence of p^, we use band structure param- 
eters obtained from ARPES and interlayer resistance data on the two-layer cuprate T^2201. The 
ARPES dataiasl] can be reasonably fit by a tight binding model with nearest and next-nearest hop- 
ping parameters: = -2?[cos + cos k^] - At' cos k^ cos ky - with k measured from {n/a, n/a), 
t' It = 0.42 and ^a/t = 1.36. The resulting Fermi surface is shown in Fig. [Tl In this material 
the interlayer hopping parameter t±{k^, ky) vanishes by symmetry at 8 points on the Fermi sur- 
face (along kx = ky and k^ = directions). It can be modeled (according to AMRO dataUJ]) as 
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FIG. 2: The non-monotonic dependence of the interlayer resistivity p„ on the magnitude of a J-wave 
pseudogap. Main panel: The vertical axis is the interlayer resistivity, averaged over the direction of the 
in-layer magnetic field (I)b, in units of the normal state zero-field value po. The horizontal axis is the 
magnitude of the J-wave pseudogap and the different curves are for different field strengths. For weak 
fields (a>cT << 1), turning on the gap has no effect other than reducing the carrier density so the resistance 
increases with KqIUbT . In strong fields, the opening of a small gap reduces the average quasiparticle velocity 
and the associated Lorentz force responsible for the large magnetoresistance. This effect results in an initial 
drop in the interlayer resistance. As Aq /ksT becomes large, the reduction of carrier density eventually 
overrides this effect and p^^ begins to increase. Upper inset: The pseudogap density of states v{E) in terms 
of the normal state value vq. The scattering rate r"' depends on E through the density of states. 

t±W = ?±[sin20 -I- ^6 sin 60 -I- (k(, - 1.0) sin 100] with = - l.O and = 0.71. The energy 
scale tj_ can be absorbed into the zero-field, normal-state resistivity but the large anisotropy in tj_{(p) 
contributes to the strong ^g-dependence observed for this material in the normal state. Moreover, 
since vanishes at the nodes, the magnitude of p^^ becomes extremely large in the nodal limit 
Ao/ keT » 1. 

The anisotropic magnetoresistance in the normal state is illustrated in the polar plots of 
Pzzk4>B)IPo versus (ps in Panel A of Fig. HJ A field of ojcT ~ 0.5 is sufficient to reveal the strong 



13 




Ajj/kgT=0, vary cOj^x 




Aj^/kgT=10, vary co^x 





co^x=0.5, vary co^x=0.5, vary A^^/k^T 



FIG. 3: Anisotropy of the interlayer resistance p^^ in tlie nomal and tZ-wave pseudogap states. Solid curves 
are polar plots of Pzz/po versus (f>B where a crystal axis is along the horizontal and the band structure of 
T£ 2201 has been used. The dashed curve is the unit circle (unseen in panel B, where the radial scale 
is much larger). Panel A: The normal state (constant Ao/ksT = 0) for varying field strength; the solid 
curves from inside out are for: ojct = 0.1,0.2,0.3,0.4,0.5. Panel B: The low-temperature pseudogap state 
(Ao/kBT = 10) for varying field strength; the solid curves from inside out are for: ojct = 0.5, 1, 2, 4, 6. Panel 
C and D: A large applied field (constant a>cT = 0.5) with varying pseudogap magnitude. The solid curves 
in C are, from outside in: Ao/ksT = 0,0.2,0.4,0.6 and in D, from inside out: Ao/ksT = 0.6, 1.0, 1.4, 1.8. 
The resistance first decreases (in C) then increases (in D) as the pseudogap grows. The angle-dependence 
is reduced by the opening of a pseudogap and is eventually replaced by that associated with the anisotropic 
gap itself. 



anisotropy of the underlying band structure. Note also that the (^g-averaged magnitude of Pzz/Po 
decreases as the scattering rate increases for a given field strength. In Panel B, the nodal limit 
Ao/ksT » 1 of Pzz/Po is depicted. Here the current is coming entirely from momenta near the 
nodes and thus provides no information about the normal state band parameters elsewhere on the 
Fermi surface. The anisotropy, which has been discussed in Ref. |29|, results from unequal, and 
0B-dependent, contributions from different nodes owing to the large ratio of v//v2 where v/ is the 
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Fermi velocity is the 'gap' velocity. 

Panels C and D of Fig. [3]describe the effect that turning on a (i-wave pseudogap has on PzzispB) 
in a relatively strong field (ojcT = 0.5). The scattering rate was evaluated in the unitary limit, using 
the rounded density of states plotted in Fig. [2l In Panel C, which shows small values of Aq /ksT, 
the magnitude of p^^ decreases as the gap opens. In Panel D, which shows larger values of Aq/ ksT, 
the 0B-averaged resistance has already reached its minimum value, depicted in Fig. [21 and is thus 
growing with AqIUbT. 

It is seen, by comparing Panels A and C, that the initial effect of a small gap on Pzz((f>B) is 
similar to the effect of an enhancement in the scattering rate. The reason for this similarity fol- 
lows from the discussion of Fig. [1} the pseudogap increases the band of k-points that contribute 
to the interlayer current just as would an increase in r"^ The manner by which the (pB depen- 
dence changes as the pseudogap continues to grow in magnitude is, however, very different from 
that resulting from an increase in the scattering rate. Not only does the magnitude of p^^ vary 
non-monotonically with Ao/ksT, but Pzz{4>b) evolves to incorporate the anisotropy of Ak along 
with that already coming from fx(0) and the intra-layer band parameters. Given the success of 
the angle-dependent interlayer resistance technique in extracting precise values for several band 
structure parameters, it appears that this technique should be equally capable of obtaining both 
the magnitude and anisotropy of a pseudogap as it emerges in the over- or near optimally doped 
cuprates. 

CONCLUSIONS 

Measurements of the interlayer resistivity in layered metals, made in a magnetic field with 
varying orientation, can be used to characterize anisotropic properties within individual layers. 
Among these properties, band structure parameters and the inelastic scattering rate in various 
systems have already been extracted by this method. In this article we have extended the analysis 
of such measurements to incorporate a pseudogap with ^^-wave symmetry. A general expression 
for the interlayer resistivity in the pseudogap state was obtained via a semiclassical calculation. 

For a field along the layers, the main effect of a small pseudogap is to smooth the dependence 
of the resistivity on the in-layer field direction (pB- This occurs because, while electrons only con- 
tribute to the normal state interlayer current if they are located at a particular point on the Fermi 
surface, so that they have a velocity parallel to the magnetic field, quasiparticles with an energy 
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smaller than the pseudogap can contribute to the interlayer current from anywhere on the Fermi 
surface. The average magnitude of the interlayer resistivity first decreases, then subsequently in- 
creases as the pseudogap opens, reaching a minimum value when the magnitude of the pseudogap 
is comparable to the temperature. We hope that this work will stimulate new experiments and anal- 
ysis to detect the presence and map the anisotropy of a pseudogap in layered strongly correlated 
materials. 

We thank Matthew French and Nigel Hussey for helpful discussions and for sharing with us 
unpublished data. This work has been supported by the Australian Research Council. 
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